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are  compared  to  those  of  Kazakh  explosions.  P-waves  to  NORSAR  from  these 
earthquakes  exhibit  very  high  frequency,  and  their  spectral  corner  frequen- 
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INTRODUCTION 


Discrimination  parameters  have  been  applied  in  detail  at  only  a few  nu- 
clear test  sites  and  earthquake  regions.  This  study  is  part  of  a series  which 
will  extend  these  discrimination  studies  to  other  regions  of  shallow  earth- 
quakes. This  particular  report  examines  eleven  earthquakes  in  the  Northern 
Pamirs  which  had  low  reported  Mg  for  their  n^,  such  that  these  events  fall 
close  to  the  explosion  population  on  a M^-m^  graph.  Most  of  the  earthquakes 
in  the  Pamirs  are  probably  associated  with  large  thrust  faults,  which  are  a 
result  of  the  Indian-Eurasian  collision.  Actual  seismogram  analysis  revealed 
that  the  events  chosen  were  too  small  to  determine  their  fault-plane  solutions. 
Paramount  in  this  study  is  the  determination  of  average  Mg  and  m^  for  selected 
events.  Further,  we  investigate  other  common  discrimination  parameters  such 
as  first  motions,  corner  frequency,  complexity,  and  short-period  spectral 
ratios.  Results  for  the  Pamir  events  will  be  compared  to  Kazakh  explosions 
and  related  to  the  available  geophysical  data  on  the  crust  and  upper  mantle 
of  the  region. 
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TKCT0N1C  SETTING 

General  Features 

Figure  1 is  a map  of  the  region  of  interest.  The  area  just  to  the  west 
of  this  map  is  denoted  the  "Garm  region"  by  the  Soviets;  it  has  the  strongest 
seismicity  in  the  U.S.S.R.  (Keilis-Borok,  I960)  and  recently  has  been  notable 
for  earthquake  prediction  efforts  (Semyenov,  1969).  This  study  will  examine 
crustal  seismicity  in  an  area  east  of  Garm  in  the  Northern  Pamirs  where  seismi- 
city is  also  high.  The  seismicity  in  this  area  is  shallow  (<  70  km)  and  runs 
parallel  to  the  E-W  trend  of  the  mountain  chain  and  is  probably  a result  of 
the  ongoing  Indian-Eurasian  collision  (Molnar  and  Tapponnier,  1975).  South 
of  the  Pamirs  in  the  Hindu  Kush,  a vertical  zone  of  intermediate-depth  seismic 
activity  suggests  that  under thrust ing  of  the  lithosphere  occurred  there  in  the 
last  10  million  years  (Molnar  et  al.,  1973).  Measurements  of  vertical  and 
horizontal  movements  suggest  that  uplift  and  thrust  faulting  are  still  con- 
tinuing in  the  Pamirs. 

The  Pamirs  consist  of  a high  plateau  of  folded  Cenozoic  sediments.  A 
large  negative  Bouguer  gravity  anomaly  of  -450  mgal  in  the  Northern  Pamirs  is 
associated  with  relatively  large  crustal  thicknesses  (55-65  km)  in  that  region 
(Kosminskaya  et  al.,  1958). 


Keilis-Borok,  V.  1.  (i960).  Investigation  of  the  Mechanism  of  Earthquakes,  Eng- 
lish translation,  American  Geophysical  Union,  Washington,  D.C. 

Semyenov.  A.  N.  (1969),  Variations  in  the  travel  time  of  transverse  and  longi- 
tudinal waves  before  violent  earthquakes,  Izvestiva,  Phys.  Solid  Earth, 

No.  4,  245. 

Molnar,  P.,  and  P.  Tapponnier  (1975),  Cenozoic  tectonics  of  Asia:  effects  of 
a continental  collision.  Science , 189,  419. 

Molnar,  P.,  T.  Fitch,  and  F.  Wu  (1973),  Fault  plane  solutions  of  shallow  earth- 
quakes and  contemporary  tectonics  in  Asia,  Earth  and  Planetary  Science 
Letters,  19,  101. 

Kosminskaya,  I.,  G.  Mikhota,  and  V.  Tulina  (1958),  Crustal  structure  of  the 
Pamir-Alai  zone  from  seismic  depth-sounding  data,  Izvestiva,  Geophysics 
Series,  673. 
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Source  Mechanisms  of  Earthquakes 

NE1S  epicenters  for  events  in  this  area  from  January  1971  through  Febru- 
ary 1975  are  shown  in  Figure  1.  The  geologic  faults  in  this  figure  were  in- 
ferred from  a satellite  photomosaic  of  the  area  and  published  geologic  or 
tectonic  maps.  Most  of  the  earthquakes  are  associated  with  the  Pamir  thrust 
fault  and  the  other  thrust  faults  which  lie  parallel  to  it  (Shirokova,  1967). 
In  the  northern  part  of  this  region,  some  of  the  events  lie  close  to  a left- 
lateral  strike-slip  fault.  The  compressive  stress  component  for  events  along 
the  Pamir  thrust  fault  is  generally  directed  perpendicular  to  the  structural 
trend  of  the  Pamir  Mountain  system  (Kitsema,  1966).  Molnar  et  al.  (1973) 
found  both  thrust  faulting  and  strike-slip  faulting  associated  with  the  Pamir 
thrust.  In  general,  however,  the  pattern  of  fault  plan  solutions  is  very  com- 
plex and  unpredictable  and  has  little  correlation  with  known  faults  (Keilis- 
Borok,  1960). 

Velocity  Model 

Seismic  studies  in  the  Pamirs  report  an  average  crustal  thickness  of 
> 50  km  (Kosminskaya  et  al.,  1958)  wLth  depths  to  65  km  indicated.  Table  1 

shows  a general  crustal  and  upper  mantle  velocity  structure  derived  from  deep 
seismic  sounding  in  the  area  of  interest;  the  actual  structure  varies  rapidly 
over  the  area  encompassing  the  earthquakes  we  have  selected  for  this  study 
though.  Seismological  data  show  that  a slight  low-velocity  zone  with  Vp  - 


Shirokova,  E.  (1967),  General  features  in  the  orientation  of  principal  stres- 
ses in  earthquake  foci  in  the  Mediterranean-Asian  seismic  belt,  Izvestiya, 
Physics  of  the  Solid  Earth,  12,  12. 

Ritsema,  A.  (1966),  The  fault-plane  solutions  of  earthquakes  of  the  Hindu- 
Kush  center,  Tectonophysics , 2^,  147. 

Molnar,  P.,  T.  Fitch,  and  F.  Wu  (1973),  Fault  plane  solutions  of  shallow  earth- 
quakes and  contemporary  tectonics  in  Asia,  Earth  and  Planetary  Science 
Letters,  19,  101. 

Keilis-Borok,  V.  I.  (1960),  Investigation  of  the  Mechanism  of  Earthquakes,  English 
translation,  American  Geophysical  Union,  Washington,  D.C. 

Kosminskaya,  I.,  G.  Mikhota,  and  Y.  Tulina  (1958),  Crustal  structure  of  the 
Pamir-Alia  zone  from  seismic  depth-sounding  data,  izvestiya  Geophysics 
Series,  673. 
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7.9  km/ sec  extended  from  10  to  20  km  below  the  Mo ho  to  about  180  km  depth 
(Aliev  et  al.,  1976,  and  Vinnik  and  Lukk,  1974).  This  contrasts  with  the  much 
more  pronounced  low-velocity  zone  under  the  western  United  States  where  Vp  - 
7.7  km/sec  (Archambeau  et  al.,  1969).  There  is  evidence  from  pP  reflections 
(Vinnik  and  Godzikovskaya,  1973)  of  a zone  of  higher  than  normal  upper-mantle 
attenuation  to  the  south  and  east  of  the  area  encompassing  the  earthquakes  of 
this  report,  but  this  same  data  indicated  that  attenuation  is  less  than  normal 
under  the  Pamirs  itself. 


Aliev,  S.,  N.  Beliaevsky,  E.  Butovskaya,  B.  Volvovsky,  I.  Volvovsky,  0.  Kras- 
nopevtseva,  V.  Pak,  M.  Polshkov,  V.  Rubai lo,  V.  Sallogub,  B.  Tal-Virsky, 

E.  Tregub,  1.  Khamrabayev,  and  G.  kharechko  (1976),  The  seismic  experi- 
ment in  the  Northern  Pamirs,  in  Geodynamics:  Progress  and  Prospects , 
ed.  G.  Drake,  American  Geophysical  Union. 

Vinnik.  L.  P.,  and  A.  A.  Lukk  (1974),  Lateral  inhomogeneities  of  the  upper  mantle 
under  the  Pamir  and  Hindu  Kush,  lzv. , Earth  Physics,  No.  1,  9. 

Archambeau,  C.,  E.  Flinn,  and  1).  Lambert  (1969),  Fine  structure  of  the  upper 
mantle,  JL  Geophys,  Res . , 74 , 5825. 

Vinnik,  L.  P.,  and  A.  A.  Godzikovskaya  (1973),  Lateral  variations  of  the  absorp- 
tion by  the  upper  mantle  beneath  Asia,  lzv. , Eartli  Physics,  No.  1,  3. 
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DATA 

Earthquake  Selection 

For  this  study  earthquakes  were  selected  which  had  low  reported  for 
their  m^  so  that  these  events  fall  close  to  the  explosion  population  on  an 
Ms-mb  KraP^*  The  NE1S  list  suggests  both  shallow  and  Intermediate  depth  events 
in  this  region.  Most  of  the  intermediate  depth  activity  takes  place  south  of 
the  area  of  interest  in  the  Hindu  Kush.  Shallow  depth  events  (less  than  70 
km)  predominate  in  the  Northern  Pamirs.  Most  of  tin  seismic  activity  in  the 
Northern  Pamirs  is  associated  with  the  Pamir  thrust  tault,  as  shown  in  Figure 
1.  The  earthquakes  chosen  were  limited  to  tin*  years  1471  to  1475  so  that  we 
could  utilize  the  data  from  the  large  seismic  arrays  and  the  HG1.P  network. 

We  selected  a total  of  11  earthquakes  as  listed  In  Table  11  and  shown  in 
Figure  1.  Table  II  also  lists  the  depths  determined  from  pP  observations  at 
LASA  and  NORSAK.  All  of  the  events  are  located  within  the  crust.  Events  5 
through  11  are  all  closely  located,  and  although  no  pP  u^ta  was  available  at 
LASA  or  NORSAK  for  event  10,  a pP  observation  at  SHI  implies  a focal  depth  of 
25  km. 

Explosion  Selection 

The  Pamir  earthquakes  were  compared  to  the  Kazakh  explosions  which  lie  to 
the  northwest  and  northeast  of  the  Pamirs.  We  selected  a total  of  10  explo- 
sions as  listed  in  Table  III. 

Seismic  Stations 

Digital  data  from  the  three  arrays  ALl’A,  LASA,  and  NORSAR  and  from  the 
available  HGLP  stations  and  film  data  from  selected  WWSSN  stations  were  ga- 
thered for  these  events.  The  LASA  AO  subarray  and  NORSAK  C3  subarray  short- 
period  P recordings  are  shown  in  Figure  2.  WWSSN  stations  were  selected  on 
the  basis  of  magnification  and  proximity  to  the  region  of  study,  and  it  is 
unlikely  that  the  stations  not  used  here  would  significantly  add  to  the  data 
base  for  these  Pamir  events. 
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LASA  and  NORSAR  short- period  P recordings  for  the  l’amir  events 
studied 
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SIONAL  ANALYSIS 


Source  effects 

We  wanted  Initially  to  Identity  a source  mechanism  for  each  event  so 
that  we  could  predict  the  radiation  pattern  for  body-wave  and  surface-wave 
phases.  Determination  of  the  fault  planes  for  earthquakes  whose  magnitudes 
are  below  b.O  is  generally  not  reliable  witli  teleselsmlc  data,  and  all  our 
earthquakes  have  m^  less  than  6.  Few  clear  first  motions  were  recorded  from 
the  data  we  had  collected.  ISC  bulletin  data  was  available  for  events  1 
through  4;  only  events  .1  and  4 had  sufficient  numbers  ol  first  motions  to 
attempt  a focal  mechanism  plot.  We  have  plotted  the  ISC  (International  Seis- 
mo logical  Center)  first  motions  for  events  3 and  4 in  Figure  3.  Unfortunately, 
the  short-period  data  were  not  consistent  and  It  is  impossible  on  the  basis  ol 
the  plots  alone  to  determine  fault  planes  because  dllatational  and  compres- 
sional  first  motions  do  not  separate. 

Corner  frequencies  and  seismic  moments  for  our  I’amir  events  have  been  es- 
timated from  short-period  1j\SA  and  NORSAR  spectra.  The  spectra  are  shown  in 
Figure  4.  These  spectra  are  from  the  phased  beams  of  the  AO  short-period  sub- 
array at  LASA  and  Cl  short-period  subarray  at  NORSAR.  The  sample  length  was 
6.4  seconds.  The  signals  have  been  tapered,  and  the  Instrument  response  was 
removed  from  the  signal  and  noise  spectra;  but  noise  spectra  have  not  been 
subtracted  from  the  signal  spectra.  Attenuation  was  removed  from  the  spectra 

by  multiplying  by  the  factor  explnft*]  with  a t*  ol  .53  for  a ui  " source  model 

-3  -2-3 

and  .34  tor  ui  source  models  at  l.ASA  and  t*  ot  .00  for  u>  and  u>  source 

models  at  NORSAR.  The  basis  for  these  t*  values  will  be  shown  later  in  this 

report.  Corner  frequencies  were  estimated  with  the  assumption  ot  complete 

-2  -3 

stress  drop  over  a circular  or  square  iault  plane  and  a io  or  io  asymptotic 
relation  at  high  frequencies.  Note  that,  for  many  of  the  spectra  ot  Figure  4, 
dashed  lines  are  used  to  indicate  the  long-period  level  and  the  corner  fre- 
quency. This  Indicates  our  lack  of  confidence  in  the  estimates.  Occasionally 
we  computed  spectra  for  25.6-second  windows  for  those  cases  where  the  slgnal- 
to-nolse  ratio  was  high.  For  these  cases  the  shape  ol  the  spectra  were  suf- 
ficiently Identical  to  those  for  the  6.4-second  windows  that  we  lujvo  not 
shown  them. 
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DISCRETE  NOISE  AMPLITUDE  SPECTRUM 
DISCRETE  SIGNAL  AMPLITUDE  SPECTRUM 


A AO  and  NORSAR  C3  subarray  spectra  of  P waves  from  Pamir 
thquakes  with  instrument  response  and  attenuation  removed 


DISCRETE  NOISE  AMPLITUDE  SPECTRUM 
DISCRETE  SI6NAL  AMPLITUDE  SPECTRUM 


'•’ORSAR  C3  subarray  spectra  of  P waves  from  Pamir 
• ith  instrument  response  and  attenuation  removed 


DISCRETE  NOISE  AMPLITUDE  SPECTRUM 
DISCRETE  SIGNAL  AMPLITUDE  SPECTRUM 


A AO  and  NORSAR  C3  subarray  spectra  of  P waves  from  Pamir 
thquakes  with  instrument  response  and  attenuation  removed 


DISCRETE  NOISE  AMPLITUDE  SPECTRUM 
DISCRETE  SIGNAL  AMPLITUDE  SPECTRUM 


FREQUENCY  (Hz) 

A AO  and  NORSAR  C3  subarray  spectra  of  P waves  from  Pamir 
thquakes  with  instrument  response  and  attenuation  removed 


DISCRETE  NOISE  AMPLITUDE  SPECTRUM 
DISCRETE  SI6NAL  AMPLITUDE  SPECTRUM 


NORSAR  C3  subarray  spectra  of  P waves  from  Pamir 
with  instrument  response  and  attenuation  removed 


DISCRETE  NOISE  A 


A AO  and  NORSAK  C3  subarray  spectra  of  P waves  from  Pamir 
thquakes  with  instrument  response  and  attenuation  removed 


Figure  4 (cont.)  LASA  AO  and  NORSAR  C3  subarray  spectra  of  P waves  from  Panir 
earthquakes  with  instrument  response  and  attenuation  removed, 


DISCRETE  NOISE  AMPLITUDE  SPECTRUM 
DISCRETE  SIGNAL  AMPLITUDE  SPECTRUM 


A AO  and  NORSAR  C3  subarray  spectra  of  P waves  from  Pamir 
thquak.es  with  instrument  response  and  attenuation  removed. 


DISCRETE  NOISE  AMPLITUDE  SPECTRUM 


NORSAR  C3  subarray  spectra  of  P waves  from  Pamir 
with  instrument  response  and  attenuation  removed 


DISCRETE  NOISE  AMPLITUDE  SPECTRUM 
DISCRETE  SI6NAL  AMPLITUDE  SPECTRUM 


FREQUENCY 

Figure  4 (cont.)  LASA  AO  and  NORSAR  C3  subarray  spectra  of  P waves  from  Pamir 
earthquakes  with  instrument  response  and  attenuation  removed 


DISCRETE  NOISE  AMPLITUDE  SPECTRUM 


A AO  and  NORSAR  C3  subarray  spectra  of  P waves  from  Pamir 
thquak.es  with  instrument  response  and  attenuation  removed. 


DISCRETE  NOISE  AMPLITUDE  SPECTRUM 
OISCRETE  SIGNAL  AMPLITUDE  SPECTRUM 


Figure  4 (cont.)  LASA  AO  and  N’ORSAR  C3  subarray  spectra  of  P waves  from  Pamir 
earthquakes  with  instrument  response  and  attenuation  removed 
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Seismic  moment  was  calculated  from  the  long-period  body-wave  spectral 
displacement  levels  |*:J  estimated  in  Figure  4 using  the  relation 

4npa^D|W  | 

M - j — 

Ro* 

where  R_  was  assumed  to  be  unity  due  to  lack  of  knowledge  of  the  focal  mech- 

' ▼*  j 

anism,  values  of  pa  from  Table  I appropriate  to  the  source  depth  were  used, 
and  the  divergence  factor  D was  applied  (Ben-Menahem  et  al.,  I9u3).  A graph 
of  moment  versus  corner  frequency  is  shown  in  Figure  5.  Our  events  fall  be- 
tween the  constant  stress  drop  lines  of  Hanks  and  Thatcher  (1972),  indicating 
that  our  events  are  of  intermediate  stress  drop. 

We  attempted  to  find  the  LR  radiation  pattern  for  each  event  by  plotting 
the  ant l logs  of  the  M values  for  each  event  as  listed  in  Table  IV.  The  re- 
sults from  the  largest  M earthquake,  event  6,  are  shown  in  Figure  6.  The 
observed  distribution  of  amplitudes  does  not  readily  conform  to  any  quadripole 
radiation  pattern,  which  is  probably  due  to  varying  propagation  effects  coup- 
led with  the  large  epicentral  distances  and  the  small  magnitudes  of  the 
events  studied.  The  scatter  of  data  points  for  all  the  events  was  so  poor 
that  no  attempt  was  made  to  find  the  best-fitting  radiation  patterns. 


Fropugat  ton  Kfl'ects 

We  measured  relative  mantle  attenuation  along  the  paths  from  the  Pamirs 
to  LASA  and  NORSAR,  the  only  two  sites  for  which  we  had  short-period  digital 
data.  Our  measurement  of  attenuation  is  the  factor  t*,  which  is  the  travel 
time  along  the  rav  path  divided  by  the  average  attenuation  factor  Q. 


We  determine  t*  by  firs:  summing  and  then  normalizing  the  observed  (t*  = 
0.0)  spectra  shown  in  Figure  4 for  the  two  paths — Pamirs  to  LASA  and  Pamirs 
to  NORSAR.  These  averages  are  shown  in  Figure  7 and  8.  If  we  assume  that  the 
source  spectra  has  a f *"  falloff  at  high  frequencies,  then 


A(f)  - S*f 


-it  ft* 

•e 


Ben-Menahem,  A.,  S.  Smith,  and  T.  Teng  (1965),  A procedure  for  source  studies 
from  spectrum  of  long-period  seismic  body-waves.  Bull.  Seism.  Soc . Am., 

55,  203. 

Hanks,  T , and  W.  Thatcher  (1972),  A graphical  representation  of  seismic  source 
parameters,  Oeophys . Res.,  _77.»  4^93. 
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Seismic  moment  versus  corner  frequency  for  Pamir  earthquakes 
from  LASA  and  NORSAK  P recordings 


Table  IV 


F 


! 

* 

f 


Magnitude  Data  for 


Pamir  Earthquakes 


Event  1 


Station 

A 

% 

M 

s 

ALP 

71.6 

3.38 

ALQ 

106.3 

-4.16 

CHG 

29.8 

-4.20 

CTA 

90.0 

-3.98 

FBK 

71.5 

-4.42 

KOD 

28.5 

4.27 

KON 

44.4 

-3.65 

LAO 

95.0 

5.05 

-4.94 

NAO 

44.0 

4.39 

3.11 

POO 

20.0 

3.97 

TLO 

57.9 

-3.23 

RINGDAL 

ESTIMATE 

4.42 

3.19 

negative  sign  indicates  no  signal  visible,  and  upper 
bound  on  is  computed  from  background  noise 

Event  2 

- 


Station 

A 

% 

M 

s 

ALP 

70.9 

-3.63 

ALQ 

105.6 

4.26 

BUL 

72.6 

4.57 

CHG 

30.3 

3.90 

CTA 

90.4 

4.29 

FBK 

70.7 

4.00 

KON 

43.8 

4.57 

NAO 

43.3 

4.23 

4.56 

NUR 

36.6 

4.47 

OGD 

94.3 

4.52 

POO 

20.8 

4.28 

-4.42 

SHI 

19.6 

3.82 

-4.22 

SHL 

20.8 

3.95 

TAB 

21.0 

4.67 

3.63 

TLO 

57.5 

3.62 

RINGDAL 

4.34 

4.04 

ESTIMATE 

N 


* 

1 1 


u 


I 


Event  3 


Station 

A 

"b 

M 

s 

ALP 

71.5 

-3.46 

ALQ 

105.9 

4.61 

BUL 

71.7 

4.26 

CHG 

30.9 

3.41 

COL 

71.6 

4.95 

CTA 

91.0 

3.86 

DAG 

52.5 

4.37 

KIP 

103.4 

4.02 

KOD 

29.1 

4.33 

LAO 

94.6 

5.25 

NAO 

43.0 

4.13 

-4.14 

OGD 

94.3 

4.17 

SHL 

21.3 

3.59 

RINGDAL 

ESTIMATE 

4.41 

3.75 

Event  4 


Station 

A 

'"b 

M 

s 

ALQ 

104.6 

4.74 

ALQ 

104.6 

4.71 

BAG 

47.3 

4.17 

BUL 

73.4 

4.70 

CTA 

90.8 

5.15 

4.10 

DAG 

51.3 

4.70 

GDH 

63.5 

5.00 

KON 

42.9 

4.58 

LAO 

93.0 

5.72 

4.44 

MAT 

49.8 

4.33 

MAT 

49.8 

4.54 

4.23 

NAO 

42.4 

4.47 

4.61 

NUR 

35.6 

3.60 

OGD 

93.2 

4.49 

SDB 

78.2 

5.22 

SHI 

20.0 

4.45 

4.05 

TAB 

20.9 

3.71 

TLO 

57.0 

4.23 

RINGDAL 

4.75 

4.34 

ESTIMATE 

1 


39' 


Event  5 


Station 

A 

"b 

M 

s 

ALl’ 

70.9 

-3.53 

Al.£ 

105.7 

-4.10 

BAG 

46.5 

5.03 

4.23 

BUL 

72.9 

5.16 

CHG 

29.  H 

3.60 

CTA 

89.8 

5.00 

-4.37 

GIL 

33.3 

4.22 

KIP 

102.3 

4.28 

KON 

44.3 

4.53 

4.35 

LAO 

94.4 

4.59 

3.69 

MAT 

49.7 

4.23 

MAT 

49.7 

4.10 

MAT 

49.7 

4.05 

NAI 

52.9 

4.54 

4.58 

NAO 

43.8 

4.77 

4.07 

NOR 

36 . 9 

4.47 

3.99 

OGD 

94.7 

-3.43 

SHI 

20.0 

4.48 

TAB 

21.5 

3.43 

TLO 

58 . 0 

-4.08 

TLO 

58.0 

3.79 

R INGOAL 

ESTIMATE 

4.73 

3.92 
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Event  6 


Station 


A 


ALP 

ALQ 

ALQ 

BAG 

BUL 

CHG 

CHG 

CTA 

EIL 

GDH 

1ST 

KIP 

KON 

KON 

LAO 

MAT 

MAT 

MAT 

NAI 

NAO 

NUR 

OGD 

SHI 

SNG 

STU 

TAB 

TLO 

TLO 


71.0 
105.8 
105.8 

46.5 
72.9 
29.8 

29.8 

89.8 

33.2 

64.9 

34.0 
102.3 

44.3 

44.3 

94.5 
49.7 
49.7 

49.7 

52.8 

43.8 

36.9 
94.7 

20.0 

40.1 

46.0 

21.4 

58.0 
58.0 


5.30 

5.31 

4.35 

5.H 


4.2 

4.7- 


4.7l 

4.7 

4.4' 

4.5- 

4.76 


-3 

4 

4 

4 

4 

4 


RINGDAL 

ESTIMATE 


4.75 


Event  7 


Station 

A 

mb 

M 

s 

ALP 

70.8 

3.57 

ALQ 

105.6 

4.65 

ALQ 

105.6 

4.10 

BAG 

46.8 

5.16 

BUL 

72.9 

4.65 

CHG 

30.1 

3.84 

EIL 

33.1 

3.70 

GDH 

64.6 

5.00 

4.44 

1ST 

33.7 

3.77 

KIP 

102.3 

-4.69 

KON 

44.0 

4.94 

4.90 

KON 

44.0 

4.20 

LAO 

94.2 

5.49 

4.07 

MAT 

49.9 

3.91 

MAT 

49.9 

-4.11 

MAT 

49.9 

4.12 

NAI 

52.8 

4.66 

NAO 

43.5 

5.08 

4.30 

NOR 

36.6 

4.80 

3.81 

OGD 

94.4 

-4.85 

SHI 

20.7 

4.24 

STU 

45.7 

4.80 

TAB 

21.2 

4.67 

3.49 

TLO 

57.7 

4.19 

TLO 

57.7 

4.08 

RINGDAL 

ESTIMATE 

4.78 

4.06 
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Event  8 


Station 

A 

% 

M 

s 

ALP 

71.0 

-3.48 

ALQ 

105.8 

-4.62 

ALQ 

105.8 

4.21 

BAG 

46.5 

5.08 

BUL 

72.9 

4.74 

CHG 

29.8 

4.26 

CHG 

29.8 

4.30 

3.96 

CHG 

29.8 

3.66 

CTA 

89.8 

4.46 

EIL 

33.2 

4.20 

KIP 

102.4 

4.27 

KON 

44.3 

4.42 

LAO 

94.5 

4.15 

MAL 

59.8 

4.39 

MAT 

49.8 

4.13 

MAT 

49.8 

4.36 

4.14 

NAO 

43.9 

4.77 

4.59 

NUR 

36.9 

4.46 

3.99 

OGD 

94.8 

3.42 

SHI 

20.0 

3.82 

-4.24 

TAB 

21.5 

4.90 

3 . 60 

TLO 

58.0 

4.56 

TLO 

58.0 

4.32 

RINGDAL 

ESTIMATE 

4.56 

4.10 
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Event  11 


< 


li 

1 

<* 


Station 

A 

"b 

M 

s 

ALP 

70.9 

3.65 

AL£> 

105.8 

-4.67 

ALQ 

105.8 

4.94 

BAG 

46.5 

4.16 

BUL 

72.9 

4.56 

CHG 

29.8 

4.09 

4.03 

CHG 

29.8 

3.82 

CTA 

89.8 

4.08 

DAG 

52.5 

4.70 

EIL 

33.3 

4.01 

1ST 

34.0 

3.83 

KIP 

102.3 

4.31 

KON 

44.2 

-5.11 

KON 

44.2 

4.72 

LAO 

94.4 

5.02 

-4.43 

LEM 

55.7 

-4.49 

MAT 

49.7 

4.38 

MAT 

49.7 

4.67 

4.41 

NAI 

52.9 

4.15 

NAO 

43.8 

4.54 

NUR 

36.9 

4.80 

-4.77 

PMG 

83.3 

5.60 

SHI 

20.0 

4.20 

4.29 

SNG 

40.1 

3.57 

TLO 

58.0 

4.11 

TLO 

58.0 

3.88 

RINGDAL 

ESTIMATE 

4.70 

4.14 
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whore  A < f ) Is  the  observed  spectrum  (corrected  for  seismograph  response),  S 
is  the  source  spectrum  scale  constant,  and  f is  frequency  in  Hz.  Then  after 
taking  the  natural  log  of  this  equation 

Ln l A ( f ) J + 2*ln(f)  - ln(S)  — lift* 

or 

ln(A)  + 2»ln(f)  — nt*f  + ln(S) 

If  we  plot  ln(A)  + 2»ln(f)  versus  frequency  for  each  path,  as  shown  in  Figures 
9 and  11,  the  slope  of  the  graph  is-nt*.  We  calculated  the  slope  in  the  fre- 
quency range  1.0  to  2.5  Hz,  where  the  signal-to-noise  ratio  was  highest,  and 

presumably  beyond  the  corner  frequency  so  that  the  high-frequency  asymptotic 
-2 

slope  of  f characterizes  the  source  spectrum.  Similarly,  the  equation  for 
the  source  spectra  with  an  f ^ falloff  is 

ln(A)  + 3*ln(f)  =-nt*f  + ln(S) 

If  we  plot  the  left-hand  side  of  this  relation  versus  frequency  for  each  path, 

as  shown  in  Figures  10  and  12,  the  slope  of  the  graph  is  again-nt*;  but  t*  now 

-3  -2 

is  based  on  an  assumed  f source  spectrum.  The  t*  values  for  an  f source 

model  are  .53  + .09  for  the  Pamir  to  LASA  path  and  .00  + .05  for  the  Pamir  to 
NORSAR  path.  The  t*  values  for  an  f source  model  are  . 34  + .08  for  the  Pam- 
ir to  LASA  path  and  -.18  + ,0b  for  the  Pamir  to  NORSAR  path.  Since  a negative 
t*  path  is  unrealistic,  the  f ^ t*  value  for  the  Pamir  to  NORSAR  path  is  set 
equal  to  zero.  The  Pamir  to  NORSAR  path  has  less  attenuation  as  expected 

since  the  path  is  mostly  through  continental  shield  areas.  Values  of  t*  for 

_ 2 

the  Kazakh  explosions  for  f source  models  were  also  calculated.  These  t* 
values  are  .16  + .10  for  the  East  Kazakh  (events  1 through  4 and  9 and  10)  to 
LASA  path,  .05  + .08  for  the  East  Kazakh  to  NORSAR  path,  .29  + .13  for  the  West 
Kazakh  (events  5 through  8)  to  LASA  path,  and  .19  + .04  for  the  West  Kazakh 
to  NORSAR  path.  As  expected  the  t*  values  for  the  paths  to  NORSAR  are  lower 
than  the  values  for  the  paths  to  1ASA  since  the  paths  to  NORSAR  are  through  con- 
tinental shield  areas  while  those  to  1ASA  pass  under  tectonically  active  pro- 
vinces . 


ZIM  (F> 


IM  (AMP)  + 3LN  (F) 


Von  Seggern  and  Sobel  (1975)  havo  presented  a method  of  predicting  20- 
second  Rayleigh-wave  amplitudes  by  geometrical  ray  tracing.  Their  work  shows 
that  lateral  changes  in  20-second  LR-phase  velocity  can  cause  major  multipath- 
ing  and  also  large  teleseismic  amplitude  variations  due  to  small-scale  refrac- 
tion effects.  Figure  11  illustrates  the  paths  of  rays  which  result  lor  earth- 
quakes 5 through  11,  which  are  grouped  near  19. IN,  71. 9E.  No  attempt  is  made 
to  quantitatively  relate  these  patterns  to  observed  LR  amplitudes  since  we 
cannot  correct  the  predicted  amplitudes  for  source  mechanism.  We  can  quali- 
tatively compare  the  pattern  of  rays  to  M residuals  calculated  from  the  magni- 
tude data  in  Table  IV.  The  M residuals  at  ALFA  are  negative,  as  can  be  pre- 
dicted from  the  ray  tracing  diagram  which  shows  diverging  rays  near  ALFA.  The 
M residuals  at  KIL  average  to  about  0,  as  would  be  expected  since  the  ray  pat- 
tern is  smooth  near  FIL.  Where  the  rays  cross,  we  can  expect  either  destruc- 
tive interference  (large  negative  residuals  as  at  CHG,  1ST,  TAB,  and  TLQ  in 
Table  IV)  or  a mixture  of  both  constructive  and  destructive  interferences  (as 
at  CTA,  LASA,  and  MAT) . Stations  which  record  a stable  amplitude  are  desirable 
for  calibrating  a source  region  for  amplitude.  For  events  from  the  Pamirs, 
such  stations  could  be  located  in  eastern  China,  southern  India,  western  Saudi 
Arabia,  Iraq,  Israel,  Finland,  and  northcentral  Russia.  The  complicated  LR 
ray-tracing  figure  partially  explains  why  we  were  unable  to  determine  LR  radi- 
ation patterns  from  our  amplitude  observations. 


von  Seggern,  I).,  and  F.  Sobel  (1975),  Experiments  in  refining  Ms  estimates  for 
seismic  events,  SDAC  Reprot  No.  TR-75-17,  Teledyne  Geotech,  Alexandria, 
Virginia. 
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DISCRIMINATION  ASPECTS 


In  this  and  the  following  sections  we  apply  several  common  discriminants 
between  earthquakes  and  explosions  to  the  Pamir  earthquakes  and  Kazakh  explo- 
sions. We  had  chosen  those  Pamir  events  which,  during  preliminary  screening, 
were  of  low  for  their  m^.  Then  we  Independently  determined  average  and 
m.  values  for  the  Pamir  events  in  Table  it  using  WWSSN  stations,  UCLP  stations 
and  the  ALPA,  LASA,  and  NOKSAK  arrays.  The  data  used  in  the  average  estimates 
are  listed  in  Table  IV.  Negative  values  in  the  table  indicate  noise  measure- 
ments which  were  used  as  upper  bounds  for  the  signal  amplitudes.  Magnitudes 
were  computed  according  to  the  formulas 

n ^ *=  logfA/T)  + B(A) 

M = log(A/T)  + 1.6b  iogA  + 0.3 

where  A = one-half  the  peak-to-peak  maximum  recorded  amplitude  reduced  to 
mp  ground  displacement, 

T = period  in  seconds  (restricted  to  17  Lo  23  sec  ior  M calculations), 

A » epicentral  distance  in  degrees, 

15(A)  “ Cutenberg-Richter  correction  term  for  surface-focus  P waves. 

We  also  independently  determined  M for  the  Kazakh  explosions  in  Table  ill  us- 
ing the  NOKSAK  array  and  HChP  stations.  The  data  used  in  the  average  estim- 
ates tor  explosion  M are  listed  In  Table  V,  and  explosion  m.  values  are  from 

s b 

NK I S list.  We  have  throughout  used  a method  ot  magnitude  estimation  proposed 
hv  Kingdal  (1*176),  in  which  the  magnitudes  at  the  individual  stations  are  as- 
sumed to  follow  a gaussian  distribution.  Among  this  distribution,  some  mag- 
nitudes will  fall  below  the  noise  level,  and  Kingdal's  method  then  substitutes 
a noise  measurement  at  those  stations  which  do  not  detect  and  computes  the  max- 
imum likelihood  estimate  ot  magnitude  based  on  measured  signals  and  noise.  The 
oiled  ot  tills  procedure  is  to  more  accurately  define  the  magnitude  of  events 
not  widely  recorded.  Based  on  this  method,  average  magnitude  values  for  small 
events,  where  many  ot  the  readings  are  noise  levels,  are  lower  than  what  would 


Kingdal,  F.  (l*»7b),  Maximum- like 1 1 hood  estimate  ot  event  magnitude,  Hu  1 1 . 
Seism,  Soc . Am. , bb,  7S‘>. 
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Magnitude  Data  for  Kazakh  Explosions 
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PAMIR  EARTHQUAKES 


thquakes  and  Kazakh 


result  if  only  measured  signal  amplitudes  were  used.  The  magnitude  averages 
for  the  Pamir  earthquakes  and  Kazakh  explosions  are  shown  in  Figure  14.  Ar- 
rows indicate  Mg  values  based  solely  on  noise  measurements.  In  these  cases 

the  maximum  M value  was  used  as  the  M value.  The  line  M *m,  - 1.5  (arbi- 
s s s b 

trarily  drawn  here)  clearly  separates  the  Pamir  earthquakes  from  the  Kazakh 
explosion  population.  This  line  would  also  separate  the  suite  of  Asian  events 
studied  by  Dahlman  et  al.  (1974).  The  Kazakh  explosions  are  located  in  two 
test  areas  approximately  10  to  15  degrees  northwest  and  northeast  of  the  Pamir 
earthquakes. 

Corner  Frequency 

Figure  15  is  a plot  of  |c^|  versus  corner  frequency  for  the  Pamir  earth- 
quakes and  Kazakh  explosions  recorded  at  LASA  and  NORSAR.  No  long-period  P 
phases  were  recorded  at  the  long-period  arrays  due  to  the  small  magnitudes  of 
the  events,  so  the  values  of  |si^|  were  estimated  from  the  spectral  data  of 
the  short-period  vertical  component.  Hanks  and  Thatcher  (1972)  have  shown 
that  for  a given  long-period  level  Aleutian  explosions  have  a higher  corner 
frequency  than  Aleutian  earthquakes.  The  physical  basis  of  this  discriminant 
is  the  smaller  source  time  dimension  of  the  explosion  for  a given  long-period 
level.  However,  Figure  15  shows  no  separation  between  the  Pamir  earthquakes 
and  the  Kazakh  explosions.  Furthermore,  there  is  only  a slight  trend  toward 
higher  corner  frequencies  for  lower  |c^|  for  both  the  explosions  and  the 
earthquakes. 

One  possible  reason  that  no  separation  was  observed  is  that  this  stu- 
dy used  only  short-period  P recordings,  whereas  Hanks  and  Thatcher  used  long- 
period  body  wave  recordings.  Another  possible  explanation  is  that  due  to  the 
low  magnitudes  of  the  events  studied  here,  a clear  asymptotic  value  of  |c^| 
was  not  possible  in  most  cases.  Finally,  the  lack  of  separation  could  be  due 
to  characteristics  of  the  sources  examined  in  this  study  which  were  not  present 


Dahlman,  0.,  H.  Israelson,  A.  Austegard,  and  G.  Hornstrom  (1974),  Definition 
and  identification  of  seismic  events  in  the  USSR,  Bull.  Seism.  Soc.  Am. , 
64,  607. 

Hanks,  T,  and  W.  Thatcher  (1972),  A graphical  representation  of  seismic  source 
parameters,  J.  Geophys.  Res. , 77,  4393. 
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in  the  Aleutian  events  studied  by  Hanks  and  Thatcher.  For  example,  the  source 
material  in  the  Aleutians  was  lava  and  bedded  tuffs  whereas  the  source  material 
for  the  Kazakh  explosions  was  presumed  to  be  granite. 


Long-Period  Body-Wave  Excitation 

There  were  too  few  long-period  P observations  to  determine  average  event 
ratios  of  these  quantities  to  long-period  LR  ground  displacement.  All  of  the 
long-period  P observations  were  close  to  the  noise  level,  as  would  be  expect- 
ed for  these  small  magnitude  events.  No  clear  long-period  S waves  were  obser- 
ved on  the  recordings  analyzed  for  the  Pamir  earthquakes. 

Depth  of  Focus 

All  of  the  Pamir  events  studied  here,  except  event  10,  had  apparent  pP 
phases  recorded  at  either  LASA  or  NORSAR  or  both  arrays.  No  LASA  or  NORSAR 
short-period  data  was  available  for  event  10;  however,  pP  for  event  10  was 
observed  at  a few  of  the  WWSSN  stations.  The  pP  phases  recorded  at  the  WWSSN 
stations  agreed  with  the  depths  determined  at  LASA  and  NORSAR  for  all  cases. 

Thus  many  low  M^-m^  events  in  this  region  could  be  identified  as  earthquakes 
on  the  basis  of  pP  observations.  pP  was,  of  course,  not  identified  for  the 
explosions  and  does  not,  therefore,  function  as  a positive  discriminant.  Al- 
so, for  crustal-depth  events  there  is  little  pP  moveout,  and  thus  positive 
identification  of  this  phase  is  not  possible  when  the  observing  stations  are 
in  a limited  epicentral  range. 

Complexities 

We  have  computed  the  "complexity"  parameter  as  seen  as  LASA  and  NORSAR 
for  the  Pamir  earthquakes  and  Kazakh  explosions  in  the  manner  given  by  Lambert 
et  al.  (1969).  Figure  16  shows  the  complexity  values  versus  m^.  LASA  beams 
for  the  Pamir  events  show  apparent  pP  signals  for  4 of  the  5 cases  where  complex- 
ity was  determined.  NORSAR  beams  show  apparent  pP  signals  for  5 of  the  8 cases 
where  complexity  was  determined.  We  believe  that  the  lack  of  a clear  pP  sig- 
nal in  the  other  cases  was  due  to  noisy  data  or  to  a pP  signal  that  could  not 
be  visually  separated  from  the  P coda.  Complexity  numbers  were  not  calculated 
for  some  cases  because  of  the  low  signal-to-noise  ratio,  which  would  have  made 
those  complexity  numbers  doubtful.  The  complexity  numbers  for  the  Kazakh  ex- 


Lambert,  D. , D.  von  Seggern,  S.  Alexander,  and  G.  Galat  (1969),  The  L0NGSH0T 
Experiment,  Volume  II.  Comprehensive  Analysis,  SDL  Report  No.  234, 
Teledyne  Geotech,  Alexandria,  Virginia. 
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Figure  16.  Complexity  versus  m,  for  Pamir  earthquakes  and  Kazakh  explo- 
sions from  LASA  and  NOKSAR  P recordings 
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plosions  were  low  (3.08-5. 6b)  as  would  be  expected  for  an  explosion.  The  com- 
plexity numbers  for  the  Pamir  events  were  generally  higher  (A. 24-11. 09) , but 
overlap  the  numbers  for  the  Kazakh  explosions.  Due  to  the  relatively  low  at- 
tenuation of  the  Pamir  earthquake  signals  at  LASA  and  NORSAR,  we  cannot  ascribe 
the  higher  complexity  of  the  earthquakes  to  attenuation  of  the  initial  P as 
suggested  by  Douglas  et  al.  (1973);  therefore  the  additional  earthquake  coda 
energy  must  be  related  to  source  mechanism  and  crustal  effects. 

Spectral  Ratios 

Spectral  ratios  have  been  calculated  at  LASA  and  NORSAR  according  to  a 
form  suggested  by  Lacoss  (1969): 

R =;1,95A(f)df  / ,*85A(f)df 
1.55  .45 

where  sums  of  the  equivalent  terms  of  the  discrete  Fourier  transform  have  re- 
placed the  amplitude  spectrum  integrals  over  A(f).  These  ratios  are  plotted 
against  m^  in  Figure  17  for  the  Pamir  earthquakes  and  Kazakh  explosions  for 

the  case  t*  = 0,  that  is,  using  the  uncorrected  spectra.  The  spectral  ratios 

_2 

for  the  same  P recordings  are  plotted  in  Figure  18  with  t*  values  for  oj 

-3 

source  models  tor  earthquakes  and  explosions  and  w source  models  for  earth- 
quakes. Bars  connect  the  cases  where  there  are  two  different  t*  values  for 
the  two  source  models.  The  amplitudes  of  the  P spectra  were  greater  than  the 
amplitudes  of  the  noise  spectra  in  the  frequency  range  0.4  to  2.5  Hz  for  all 
of  the  LASA  and  NORSAR  beams.  For  many  cases  the  P spectra  were  greater  than 
the  noise  spectra  in  the  range  0.5  to  5.0  Hz.  NORSAR  signal  beams  (C3  sub- 
array) possessed  more  high  frequencies  and  therefore  generally  had  higher 
spectral  ratios  than  the  LASA  beams  (A0  subarray),  in  accord  with  our  previous 

t*  estimates  for  the  two  paths.  The  spectral  ratio  did  not  depend  on  m,  , M , 

b s 

or  depth.  In  general,  the  Kazakh  explosions  show  larger  spectral  ratios  than 
the  Pamir  earthquakes,  but  the  explosion  and  earthquake  populations  cannot  be 
said  to  separate. 

Douglas,  A.,  P.  Marshall,  P.  Gibbs,  J.  Young,  and  C.  Blarney  (1973),  P signal 
complexity  reexamined,  Geophys,  J . , 33,  195. 

Lacoss,  R.  (1969),  A large-population  LASA  discrimination  experiment.  Tech- 
nical Note  1969-24,  Lincoln  Laboratory,  Lexington,  Massachusetts. 
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SPECTRAL  RATIO 


• PAMIR  EARTHQUAKES 
O KAZAKH  EXPLOSIONS 


mb 

Figure  17.  Short-period  P spectral  ratio  versus  for  Pamir  earthquakes 
and  Kazakh  explosions  recorded  at  hASA  and  NORSAR  for  t*  = 0 


SPECTRAL  RATIO 


"b 

Figure  18.  Short-period  P spectral  ratio  versus  for  Pamir  earthquakes 
and  Kazakh  explosions  recorded  at  LASA  and  NORSAR.  The  t* 
values^are  for  w source  models  for  earthquakes  and  explosions 
and  w source  models  for  earthquakes 
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Radiation  Pattern 


Body  waves  for  eight  of  the  eleven  Pamir  events  show  both  apparent  com- 
pressions! and  dilatational  first  motions,  suggesting  that  these  events  are 
earthquakes.  Due  to  the  low  magnitudes  of  these  events,  first  motions  were 
difficult  to  determine  on  most  records.  The  overlapping  congressional  and 
dilatational  first  motions  on  the  first  motion  diagrams  indicate  that  this  is 
not  a reliable  discriminant.  For  the  three  events  which  do  not  show  both  com- 
presslonal  and  dilatational  first  motions,  no  clear  first  motions  were  observed 
on  any  of  the  records.  The  first  motions  recorded  for  the  other  events  were 
too  few  in  number  to  determine  fault  plane  solutions  on  a focal  sphere.  We 
were  not  able  to  determine  a radiation  pattern  from  the  LK  amplitudes  owing  to 
the  large  scatter,  which  is  caused  by  propagation  effects,  such  as  1}  differ- 
ences, dispersion  effects  on  time-domain  amplitudes,  and  focusing  or  defocus- 
ing  due  to  refraction  as  illustrated  in  the  predicted  raypaths  in  Figure  13. 

S/P  Kxcltation 

As  would  be  expected  for  these  small  magnitude  events,  there  were  no 
short-period  S observations  from  the  Pamir  events  to  determine  the  ratio  of 
S/1’  ground  displacement. 
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SUMMARY 


Eleven  events  with  low  reported  for  their  ni^  from  the  Northern  Pamirs 
and  ten  Kazakh  explosions  were  examined  in  a seismic  discrimination  context. 

All  the  Pamir  events  lie  close  to  major  faults  observed  in  satellite  photo- 
graphs. It  was  impossible  on  the  basis  of  first  motions  to  determine  fault 
plane  solutions  due  to  the  small  magnitudes  of  the  events.  We  were  also  un- 
able to  determine  a radiation  pattern  from  the  l.K  amplitudes,  probably  due  to 
varying  propagation  effects.  The  following  characteristics  of  the  l’amir 
events  indicate  that  all  the  events  are  earthquakes: 

o All  the  presumed  earthquakes  had  consistent  pi’  phases 

recorded  at  LASA,  NORSAK,  or  some  of  the  WWSSN  stations. 

Although  little  moveout  for  pF  can  be  detected  for  shallow 
earthquakes,  the  clarity  and  size  of  the  phase  was  sufficient 
in  most  cases  to  identify  it  with  some  confidence. 

o All  the  Pamir  events  studied  here  fall  above  the  line 

M = nr  - 1.5.  The  Kazakh  explosions  clearly  lie  below 
s b 

the  line  M = ni^  - 1.5. 

There  are  other  areas  in  Asia,  notably  near  30°N,  95°E,  where  earthquakes 

do  not  clearly  separate  from  explosions  in  Asia  on  an  M -m^  plot.  Also  shot 

arrays  could  be  used  to  raise  the  M of  an  explosion  bv  about  0.3  M units  or 

s s 

more  without  raising  , placing  it  close  to  the  earthquake  population  studied 
here. 

Complexity  and  spectral  ratio  were  not  useful  discriminants  in  this  data 
set.  The  low  magnitudes  of  the  earthquakes  studied  here  (m^  from  4.3  to  4.8 

and  M from  3.2  to  4.4)  made  clear  discrimination  difficult  with  the  avail- 

s 

able  data  and  probably  represents  nearly  the  lowest  threshold  of  discrimina- 
tion for  the  stations  and  area  studied  in  this  report.  The  installation  of 
high-quality  SRO  stations  in  Asia  will  permit  the  examination  of  lower-magni- 
tude events  with  multistation  short-period  data. 
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